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Abstract 

The ALEPH data on the vector and axial- vector spectral functions, extracted from 
tau-lepton decays, is used in order to test local and global duality, as well as a set 
of four QCD chiral sum rules. These are the Das-Mathur-Okubo sum rule, the first 
and second Weinberg sum rules, and a relation for the electromagnetic pion mass 
difference. We find these sum rules to be poorly saturated, even when the upper 
limit in the dispersion integrals is as high as 3 GeV 2 . Since perturbative QCD, plus 
condensates, is expected to be valid for \q 2 \ > 0(1 GeV 2 ) in the whole complex 
energy plane, except in the vicinity of the right hand cut, we propose a modified set 
of sum rules with weight factors that vanish at the end of the integration range on 
the real axis. These sum rules are found to be precociously saturated by the data to 
a remarkable extent. As a byproduct, we extract for the low energy renormalization 
constant L\q the value —4L\q = 2 A3 x 10 -2 , to be compared with the standard 
value -4L 10 = (2.73 ± 0.12) x 10~ 2 . This in turn leads to a pion polarizability 
a E = 3.7 x 10~ 4 fm 3 . 
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There is a set of sum rules involving the difference between vector and axial- vector spectral 
functions, first discovered in the framework of current algebra [|lj, which are now under- 
stood as consequences of the underlying chiral symmetry of QCD. These spectral functions 
are related to the discontinuities in the complex energy plane of the two-point functions 
involving the vector and axial-vector currents 

lC(<? 2 ) = * / d'xe^ <0\T(V,(x) K t (0))|0> 

= {~9^ q 2 + qtfu) n y (g 2 ) , (1) 

n£V)=« / d'xe^ <0|T(A M (x) 4(0))|0> 

= (-Qvv Q 2 + qtfu) H-A(q 2 ) - q^qv n (g 2 ) , (2) 

where V^(x) =: q{x)^f ,j,q(x) :, A^(x) =: q{x) n f lx r y^q{x) :, and q = (u,d). In this note we will 
study specifically the chiral correlator IV_a = Uy — U^. This correlator vanishes identi- 
cally in the chiral limit (m q = 0), to all orders in QCD perturbation theory. Renormalon 
ambiguities are thus avoided. Non-perturbative contributions to this two-point function 
start at dimension d — 6, and involve the four-quark condensate. In the factorization 
approximation, and at high momentum transfer, this contribution is given by 
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The low energy behaviour of Uv~A{q 2 ) is governed by chiral perturbation theory; to one 
loop order this is given by 
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where fi n is the pion mass, f n the pion decay constant, f n = 92.4 ± 0.26 MeV 
(1 — 4/i^/g 2 ) 1 / 2 , and L w is the scale independent part of the coupling constant of the 
relevant operator in the 0(E A ) Lagrangian of chiral perturbation theory ||. The latter is 
related to the following combination of hadronic parameters 
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where < r 2 > is the electromagnetic mean squared radius of the pion, < r 2 >= 0.439 ± 
0.008 fm 2 [|j], and Fa is the axial- vector coupling measured in radiative pion decay, Fa = 
0.0058 ±0.0008 0. In the sixties, a number of sum rules were derived for the discontinuity 
of Uv-a{(1 2 ), i n the complex energy plane, from current algebra and an assumed asymptotic 
behaviour. The first sum rule is that of Das-Mathur-Okubo (DMO) || 
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W = / - [p v (s) - p A (s)) = -ft <rl> ~F A = -4 L 10 = (2.73 ± 0.12) x 10~ 2 , 
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where Pv.a{s) are related to the imaginary parts of Hy,A , and normalized according to 

M«) = ^2 t 1 + ° M) • (7) 

The second and third relations we consider are the first and second Weinberg sum rules 
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The final sum rule is a relation for the electromagnetic pion mass difference 

W 3 = J o dss In \ ±j [p v (s) - pa(s)\ = ^ (pl ± - pl«) . (10) 

Notice that this sume rule is actually independent of the arbitrary scale p, by virtue of 
the second Weinberg sum rule, Eq.(9). The above set of sum rules are exact relations in 
QCD, in the chiral SU(2) l x SU (2) r limit, i.e. for vanishing up- and down-quark masses || . 



With the advent of the ARGUS data on semi-leptonic tau-lepton decays 0, a recon- 
struction of the vector and axial-vector spectral functions became possible, albeit in the 
restricted kinematical range limited by the tau-lepton mass, thus allowing for the above 



sum rules to be confronted with experiment. This was first attempted in |I0 |. Shortly 
after, an improved parametrization of the tau-lepton decay data was discussed in [jnj , and 
used in |12| to check the first three sum rules above. Since kinematically s < M 2 < oo, 



the question is how well does one actually verify these sum rules. It turns out that by 



invoking quark-hadron duality, all three sum rules do effectively become Finite Energy 
Sum Rules (FESR). In other words, the upper limit of integration may be taken at a finite 
energy squared s = s , where the asymptotic freedom thershold s signals the end of the 
resonance region and the start of the hadronic continuum, which is well approximated by 
perturbative QCD. The conclusions of JTO j and [0 were that the ARGUS data appeared 
to saturate the sum rules reasonably well. In this paper we make use of the recent, and 
much more accurate, ALEPH data flnj on the vector and axial-vector spectral functions, 
extracted from tau-lepton decays, to test the validity of (a) local and global duality, and 
(b) the four chiral sum rules Eqs.(6),(8)-(10)(for earlier related work see e.g. [pf|). 



The concept of duality is well defined in QCD. The starting point is Wilson's Operator 
Product Expansion (OPE), extended beyond perturbative QCD by the addition of power 
corrections involving quark and gluon vacuum condensates jOJ. This OPE is expected to 
be valid in the whole complex energy plane, except in the resonance region. In the space- 
like region (q 2 < 0) the relevant scale of the OPE is determined by Aq C d — 300 — 400 
MeV, on the perturbative side, and by the scale of the condensates, also roughly 300 MeV, 
on the non-perturbative side. It is therefore reasonable to use the OPE in the region 
—q 2 > 1 GeV 2 . However, in the time-like region the relevant scale is much higher, as it is 
set by the masses of resonances in the range 1 — 2 GeV. In this region, the use of the OPE 
may be further impaired by instanton and threshold effects. Hence, the OPE may be used 
in the time-like region only at much higher scales than in the space-like region. To account 
for this fact, one distinguishes different types of duality, to wit. (i) Local duality 

Imn(s) ~ Imn(s)| QC D, (11) 

valid in QCD at very high — q 2 . (ii) Global Duality, which is an integral statement based 
on Cauchy's theorem 

P - f(s) Im U(s) ds = ~ I f(s) U(s) ds 
J Sth 7r 2m J 

- f f(s)U(s)\ QCD ds. (12) 



Alternatively, since Cauchy's theorem is valid in perturbative QCD, 
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In the above equations, is the physical threshold, sqcd the QCD threshold, and so > Sth 
is the so-called asymptotic freedom threshold. The latter must be chosen large enough so 
that II(s) ~ n(s)Q(7£) on the circle \s\ = s . The two equations above hold only if the 
weight function f(s) is analytic in the domain \s\ < s ; for instance, if f(s) is a polynomial. 
Global duality, therefore, relates weighted integrals of spectral functions (connected to 
cross sections through the optical theorem) to integrals of the QCD correlator over a large 
circle. Global duality in the interval s t h — s implies that local duality should be valid for 
s > sq. This may be seen by invoking Cauchy's theorem for a contour made of two circles 
of radii s an d Si > s , together with the connecting lines along the cut. In practice, global 
duality may often hold to a much better degree than local duality, because oscillations can 
lead to cancellations, (iii) Finally, we wish to introduce a form of global duality, let us call 
it restricted global duality, expected to set in much sooner than ordinary global duality. 
This can be accomplished if the weight functions in Eq.(12) are constrained to vanish at 
s = s . A well known example is the total hadronic decay rate of the r-lepton. Below, we 
shall make use of this form of duality to write four modified chiral sum rules, as alternatives 



A few comments on the above chiral sum rules are in order. The left hand side of Equations 
(8) and (9) can be related in a straightforward manner to contour integrals of the QCD 
correlator. The same is true of Eq.(6), provided one takes into account that 11(0) is known, 
i.e. 



In principle, duality arguments cannot be applied to Eq.(10) because of the presence of 
the logarithm. As this is only mildly singular, and it can be approximated rather well in 
the integration interval by an analytic function, one could still invoke duality. However, 
one should keep in mind that in this case duality may only be approximately valid. From 
a theoretical point of view, the four chiral sum rules are expected to be saturated for rela- 



to Eqs. (6), (8)-(10). 
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tively low values of so, on account of the very rapid fall-off of the correlator, cf. Eq.(3). 



In Fig. 1 we show the ALEPH experimental data points Jl3| on the vector and axial- 
vector spectral functions, together with our fit to these data (solid lines) and the pertur- 
bative QCD prediction (dash line). The normalization of these spectral functions has been 
changed (Eq.(7) has been multiplied by 4tt 2 ) in order to simplify comparisons with Ref. 
||13|| . Experimental errors are small, of the size of the dots in the figure, except near the 
end point. Even after taking these errors into account, the agreement between data and 
theory, which measures the validity of local duality, is not quite satisfactory. Furthermore, 
because of cancellations, the errors on the difference py — p A should be smaller than the 
vector and axial-vector spectral function errors taken in quadrature. In view of this, one 
would expect a similar unsatisfactory saturation of the four chiral sum rules, Eqs.(6), (8)- 
(10), and in fact, this is what we find. Before illustrating our results, let us write down a 
set of four chiral sum rules modified according to the principle of restricted global duality 
introduced above 
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W 3 = J Q dss In (-J [p v (s) - PA (s)} = ^ (/4 ± - /4o) , (18) 

where Eq.(15) is a combination of the DMO sum rule and the first Weinberg sum rule, 
Eq.(17) is a combination of the DMO sum rule and the first and second Weinberg sum 
rules, and the arbitrary scale in Eq.(18) has been fixed to So, which becomes the upper 
limit of integration in the four sum rules. We now show that these modified sum rules are 
saturated far better than the original sum rules Eqs.(6), (8)-(10). In Fig.2 we plot the left 
hand side (l.h.s.) of Eq.(6) computed using the fit to the data (curve(a)), and the right 
hand side (r.h.s.) (curve(b)). Agreement with the data can be considerably improved by 
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rescaling the r.h.s. of Eq.(6) from the value 2.73 x 10~ 2 to 2.43 x 10~ 2 . Figure 3 shows 
the l.h.s. of the modified DMO sum rule Eq.(15) (curve(a)) compared to the r.h.s. (curve 
(b)) after performing the above rescaling. This value of Wq has implications for the pion 
polarizability, i.e. 



where the numerical value above is the result of using the rescaled value Wq = 2.43 x 10~ 2 . 
In Fig. 4 we plot the l.h.s. of the first Weinberg sum rule Eq.(8) (curve(a)), its modified 
version, Eq.(16) (curve(b)), and their r.h.s. (curve(c)). Figure 5 shows the l.h.s. of the 
second Weinberg sum rule Eq.(9) (curve(a)), compared to its r.h.s. (curve(b)), and Fig. 
6 the corresponding curves for the modified sum rule Eq.(17). Finally, in Fig. 7 we plot 
the l.h.s. of the sum rule Eq.(lO) (curve(a)), the l.h.s. of the modified sum rule, Eq.(18) 
(curve (b)), and their r.h.s. (curve(c)). An inspection of these results clearly indicates 
that, in line with the unsatisfactory check of local duality, the original chiral sum rules do 
not appear well saturated by the data. While an overall constant rescaling of the exper- 
imental data can result in a better saturation of the DMO sum rule, this would not help 
with the other three sum rules. Hence, the problem cannot be blamed on a systematic 
overall normalization uncertainty in the data. On the other hand, by using restricted global 
duality, the four modified chiral sum rules Eqs.(15)-(18) are extremely well saturated by 
the data. 

In this note we have addressed the issue of the range of applicability of perturbative QCD, 
augmented by condensate terms (or higher twists) using the chiral correlator Ily.^ as an 
example. We have argued that the notions of duality or precocius scaling apply, in the 
space-like region, at remarkably low momentum transfers, but they must be used with great 
care when time-like momenta are involved (for discussions on the possible breakdown of 
local duality in QCD see e.g. [0-[0). In fact, even at a q 2 as large as q 2 ~ 3 GeV 2 
neither the experimental spectral function itself, nor its moments are satisfactorily given 
by their QCD counterpart. However, by taking linear combinations of moments, such that 
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they vanish at the upper limit of integration (i.e. at s = sq, the continuum threshold), a re- 
markable agreement with the QCD prediction is achieved. This we referred to as restricted 
duality. Clearly, if the (finite range) integral over the spectral function is equated to an 
integral over a large circle in the complex energy plane, the choice of polynomial integral 
kernels that vanish at the upper limit of integration causes the contribution in the vicin- 
ity of the cut to vanish. We therefore conclude that perturbative QCD may be expected 
to hold in the complex energy plane for \q 2 \ > s 0) with s as l° w as 1 GeV 2 , except near 
the cut, where the corresponding radius appears to be definitely much larger than 3.5 GeV 2 . 
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Figure Captions 



Figure 1. Our fits to the vector (curve(a)) and axial-vector (curve(b)) spectral func- 



tions together with the ALEPH experimental data JT3[. Normalization is according to: 
p VjA (s) = | [1 + {a,)] . 

Figure 2. The left hand side of Eq.(6) computed using the fits to the data (curve(a)), and 
its right hand side (curve(b)) 

Figure 3. The left hand side of Eq.(15) computed using the fits to the data (curve(a)), and 
its right hand side (curve(b)) with — 4Li rescaled to 2.43 x 10~ 2 . 

Figure 4. The left hand sides of Eq.(8) and Eq.(16) computed using the fits to the data 
(curves(a) and (b), respectively), and their right hand side (curve(c)). 
Figure 5. The left hand side of Eq.(9) computed using the fits to the data (curve(a)), and 
its right hand side (curve(b)) 
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Figure 6. The left hand side of Eq.(17) computed using the fits to the data (curve(a)), and 
its right hand side (curve(b)) 

Figure 7. The left hand sides of Eq.(lO) and Eq.(18) computed using the fits to the data 
(curves(a) and (b), respectively), and their right hand side (curve(c)). 
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